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A topological theory of the diabolical points (degeneracies) of quantum magnets is presented. 
Diabolical points are characterized by their diabolicity index, for which topological sum rules are 
derived. The paradox of the the missing diabolical points for Fes molecular magnets is clarified. A 
new method is also developed to provide a simple interpretation, in terms of destructive interferences 
due to the Berry phase, of the complete set of diabolical points found in biaxial systems such as Feg. 
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The energy levels of a quantum mechanical system 
generally tend to repel each other, so that degeneracies 
constitute exceptional events (in a sense to be specified 
below) Q. Such degeneracies, called diabolical points 
0, have recently attracted great attention in molecu- 
lar magnets 01, where they occur as a result of destruc- 
tive interference (due to the geometric Berry phase 0) 
between different tunnelling paths 0, and give rise to 
oscillations of the tunnel splitting of the ground state 
of quantum magnets 6] . Experiments on Fes molecular 
magnets have not only confirmed this 0, but also re- 
vealed the existence of further series of diabolical points, 
which, so far, could not be understood in terms of de- 
structive interference due to the geometric phase. Fur- 
thermore, some expected diabolical are missing, due to 
higher order anisotropies 0]- In this paper, a general 
topological theory of the diabolical points of quantum 
magnets is presented. Diabolical points arc character- 
ized by their diabolicity index, for which topological sum 
rules are derived. The paradox of the the missing dia- 
bolical points for Feg molecular magnets is clarified. A 
new method is also developed to provide a simple in- 
terpretation, in terms of destructive interferences due to 
the Berry phase, of the complete set of diabolical points 
found in biaxial systems such as Feg. 

The question of diabolical points, to be addressed 
in the present paper, gora back to the famous von 
Neumann- Wigner theorem 1] stating that, in a family of 
parameter-dependent hermitian Hamiltonians, acciden- 
tal degeneracies of two successive eigenvalues are found 
on submanifolds of codimension 3 of the parameter mani- 
fold. In other words, if a hermitian Hamiltonian depends 
on 3 external real parameters, such as the 3 components 
of the magnetic field, degeneracies can be found only for 
isolated values of the magnetic field, and therefore con- 
stitute a set of measure zero. Because the double-cone 
shape of the eigenenergy surfaces near such degeneracies 
resemble the toy called diaholo, they have been dubbed 
diabolical points 0. 

Interest in diabolical points has been renewed as Berry 
pointed out that they behave as magnetic monopoles 
in parameter space, i.e., that a system that is adiabat- 
ically transported around a closed circuit in parameter 
space near a diabolical point acquires a phase shift (the 



Berry phase) proportional to the solid angle of the circuit 
as seen from the diabolical point. For quantum spin sys- 
tems, it has been pointed out ^5] that the occurrence of a 
diabolical point implied by Kramers' theorem j8||, namely 
the absence of tunnelling between degenerate ground 
states of anisotropic quantum magnets of half-integer 
spin in zero field, can be understood as due to destructive 
interference between equivalent tunnelling paths whose 
Berry phase differ by an odd multiple of tt. Further, Garg 
has pointed out that, as a magnetic field is applied 
along a hard axis, the solid angle Vl enclosed between 
the two equivalent tunnelling paths joining the classical 
ground states A and B (red arrows in Fig. Ql) decreases 
from 27r to with increasing magnetic field, giving rise 
to 2J equidistant diabolical point located on the hard 
axis (red dots at Hz = in Fig. ^) . This prediction has 
been confirmed in a beautiful experiment by Wernsdor- 
fer and Sessoli ^ , who observed this oscillatory behavior 
(with 4 diabolical points on the positive hard axis) for the 
spin- 10 molecular mag net [(tacn)6Feg02(OH)i2]^+ (usu- 
ally abbreviated as Feg). Furthermore, Wernsdorfer and 
Sessoli discovered, for non-zero values of the easy-axis 
field Hz (Fig. QJi) , further series of unexpected diabolical 
points, displaying a characteristic parity alternation (red 
vs. blue points, on Fig. ^). 

Following this discovery, the complete set of diaboli- 
cal points has been identified by semi-classical, pertur- 
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FIG. 1: a,b: Schematic level diagram of a biaxial spin system 
with J = 3 and < D < is" for i?^ = (a) and H^>Q (b); c: 
diabolical points for a spin J = 3 with biaxial anisotropy; d: 
sketch of the various tunnelling paths between states A and 
B. 
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bative, or algebraic methods 0, 0, Writing the 

biaxial Hamiltonian as 7i = Tio — H • J, with Hq = 

-KJl + oi^Jl- jfj, a.nd < D < K, the diabolical 

points, corresponding to degeneracies between the states 
M and —M' (labelling as in Fig.^,b), are exactly given 

by 113 

(M - M')H°, 

- n \ i/o, (lb) 



= 
— 



(la) 



with n = 0, 1, . . . , {M + M'- 1), H° = 2y/2D{K+D) and 
= \/K^-D^. The full set of diabolical points (for 
J = 3) is shown in Fig.^. One should note that several 
diabolical points may coincide; the number of such coinci- 
dent diabolical points is the same on a given diamond, as 
indicated on Fig.^ jllj. In spite of the striking apparent 
similarity between the sets of diabolical points found on 
and off the hard axis, respectively, the latter could not be 
interpreted in terms of destructive interferences between 
tunnelling paths, which is very unsatisfactory. Further- 
more, only 4 diabolical points were observed on the posi- 
tive hard axis, instead of the 10 predicted, which has been 
explained as due to a small tetragonal anisotropy term 
[jj. However, what happens with the missing diabolical 
points remains mysterious. 

The general problem of diabolical points of quantum 
magnets may be formulated as follows. We consider a 
spin- J system with Hamiltonian 7i = Tip (J)— H • J, where 
the zero-field Hamiltonian 7io(J) is an arbitrary even 
function of the vector spin operator J. Note that the 
above Hamiltonian encompasses also the case of an ar- 
bitrary tensorial (^-factor, which can be accounted for by 
properly rescaling the field components along the prin- 
cipal axes of the ^-tensor. To discuss the properties of 
diabolical points, we take the convention to identify the 
2 J -I- 1 eigenstates by the label /i running from + J for the 
state of lowest energy to — J for the state of highest en- 
ergy, with increments of 1 (this labelling corresponds to 
the quantum number M of for H in the +z direction, 
and Ho 0). We call a diabolical point of order g (or a 
5-diabolical point), a point in H-space where g successive 
eigenstates are degenerate; the various diabolical points 
are labelled by an index z, running from 1 to Nd, the to- 
tal number of diabolical points corresponding to a given 
Hamiltonian Hq. Note that several diabolical points in- 
volving different sets of levels may coincide in H-space; 
this occurs, for instance, for the diabolical points of the 
biaxial case mentioned above. We also note, in passing, 
that the argument used by von Neumann and Wigner 
to discuss the occurrence of 2-diabolical points can be 
immediately generalized to show that the submanifolds 
on which we find the coincidence of n diabolical points 
of respective orders gi {1 < i < n) are of codimension 
d = X]"=i(5? — 1). A diabolical point where the levels /i 
to fi' (with fi' < fj,) are degenerate will be noted H^'^^^^ . 



As discovered by Berry , a quantum system in a non- 
degenerate eigenstate /x adiabatically transported around 
a closed curve C in parameter space (the external mag- 
netic field) acquires a geometric phase given by the flux 
through a surface S subtended by the circuit C of the 



Berry curvature B( 



-, where 



the sum is restricted to /i' ^ /i. The Berry curvature 
Bj-^) is divergenceless, except at diabolical points involv- 
ing the level /x, where monopole sources are located 0- 
To each diabolical point, we can associate a closed sur- 
face Si surrounding it, such that (except for coinciding 
diabolical points) no other diabolical point is enclosed 



inside E^. For a diabolical point H 



the definite- 



ness of the wavefunction implies that the flux through 
Ei of B(^) is topologically quantized, i.e., for /i2 <Ai<Mii 
Qi(^) = Jj^ B(^)-dS e Z; (by convention, we define 
Qi(^) =0 for /i < //2 or > /ii). The topological charge 
Qi(^) is known as a Chern number, an analogous to the 
Euler index of a surface in differential geometry. For a 
given diabolical point i, one can prove the following sum 
rule (minor extension of a result of 0]): '^^Qi{p) — 0. 
Further, considering a surface S enclosing all the dia- 
bolical points (one can show easily that such a surface 
exists), the flux of the Berry curvature through S, being 
a topological invariant, should remained unchanged as 
TLo is scaled down to zero, and is therefore given by the 
Chern number of a spin in a Zeeman field, which yields 
another sum rule: Qi(^) = 2/i. 

The diabolicity index of a diabolical point i for a pair 
of successive levels (/z,/i— 1) is defined as the sum of 
the topological charges up to level fx, i.e.: T^\(^)^'' = 
'^^ly^ Qi(p) - For notation convenience, for a f/-diabolical 



point with g > 2, H 



we lump the correspond- 



ing diabolicity indices into the multiplet V 



'V 



(M-l).p(M-'2) 



i (m) 



From the latter sum 



rule for the topological charges, above, we obtain the fol- 
lowing sum rules for the diabolicity indices: 



V 



(m) 



V 



E^1(;)'^ = (^ + A^)(^-(a^-1))' (2a) 



2J(J+ 1)(2J+ 1) 



(2b) 



In addition to these sum rules, the set of diabolical points, 
for a given Hamiltonian Tio , must possess all symmetries 
of Tio; in particular, the time-reversal invariance of Tio 
implies the inversion symmetry of the set of diabolical 
points. One observes immediately that all these rules 
are obeyed for the diabolical points of the biaxial sys- 
tem (Eqs. (jlal b)). where only 2-diabolical points with 
diabolicity index 1 are found. Finally, when scaling the 
Hamiltonian as TLq XHq with A > 0, the diaboli- 



cal points scale as H 



i(fj.} 



AHj|'^|, the diabohcity in- 
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FIG. 2: Spectrum and diabolical points of a spin J = 2 (a) 
and J = 5/2 (b) with cubic anisotropy: 'Hi} = EQ + K{St + Sy + 
St)/Q- Positive (resp. negative) values of the field correspond 
to a field parallel to a fourfold (resp. threefold) symmetry 
axis. The diabolical points are indicates by the solid dots. 
The corresponding diabolicity indices are indicated. 



dices remaining unchanged, whereas under reversing the 
sign of the Haniiltonian, i.e. Tip —'Ho, the diabolical 



f{p') 



H 



points and diaboUcity indices change as H^.j-^^ — ^ 

To illustrate these rules for a 



and P^r,'^ ^ -V^'^^ 



case where higher-order diabolical points occur, I show 
in Fig. 121 the spectrum and diabolical points for spins 
2 and 5/2 with cubic anisotropy, where a rich variety 
of diabolical points is obtained. The diabolicity indices 
characterize the energy dispersion near a diabolical point: 
the double-cone (diabolo) shape is obtained only for a 2- 
diabolical point of diabolicity index ±1; otherwise, a dif- 
ferent dispersion law is obtained, as seen in Fig|21 The di- 
abolicity index also influences the nature of the Landau- 
Zener tunnelling taking place at a diabolical point. A 
systematic study of these issues will be given elsewhere. 

I then address the above mentioned paradox of the 
missing diabolical points for Fes- As already indicated, 
it has been found that that the experimental observations 
are well explained quantitatively by adding to the main 
biaxial Hamiltonian a (very small) fourth order tetrago- 
nal anisotropy term H' = c(j+ -h Jl), with C < 0. 
It has been suggested that the additional anisotropy 



term might lead to a singular behavior; as we shall see, 
this explanation is both correct and incomplete. It is in- 
complete, because the sum rule (1^ would be violated 
if the diabolical would have simply disappeared. Let us 
qualitatively discuss what happens as one continuously 
switches from a biaxial anisotropy {D > 0, C = 0) to a 
tetragonal one (D = 0, C<0) (Fig. El . The effect of the 
additional term (with C < 0) is to introduce a new tun- 
nelling path (yellow arrow in Fig. ^1); for small values 
of |C| and Hx, the corresponding amplitude is negligible 
and the effect is only to displace the diabolical points 
along the hard axis, reducing the distance between the 
last ones. At a critical value of C, beyond which the am- 
plitude of the yellow path becomes larger than the one of 
the red ones, the last 2 diabolical points collide and a bi- 
furcation takes place. Beyond that point, the 2 diabolical 
points symmetrically diverge away from the x axis, to- 
wards the x-j-y axis (hard axis for tetragonal anisotropy 
with C < 0) ; this process then repeats until all diabolical 
points have moved to the x-j-y axis, for D = 0. The sce- 
nario corresponding to C > is also shown in Fig. O For 
the case of Fes (^=10), this interpretation implies that 3 
diabolical points should be located on each branch of the 
fork seen for I? > and C < in Fig.O The experimen- 
tal check of this prediction would allow to confirm the 
present topological theory of diabolical points. Finally, 
I propose the following conjecture: a spin Hamiltonian 
Ti.0 is completely determined by the set of its diabolical 
points and diabolicity indices (together with the value of 
its trace). 

I now come to the last point of this paper, namely 
the Berry phase interpretation of the diabolical points 
found at non-zero values of for the biaxial system. 
The difficulty lies in the fact that the initial (Af) and 
final state (— M') of the tunnelling paths do not gen- 
erally belong to the same set of coherent states, which 
makes the path integral approach of Ref. impracti- 




FIG. 3: Schematic representation of the evolution of diaboli- 
cal point distribution (for the two lowest states, J = 4), as the 
anisotropy progressively changes from biaxial to quadratic. 
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cable. A solution to this difficulty consists in enlarging 
the Hilbert space to comprise all possible states of a sys- 
tem of 2J spins 1/2: ji, j2, • ■ • ,j2j- The Hamiltonian 7i 
operates in this new Hilbert space by interpreting J as 
J = Y^i=i ii ■ We can ensure that the physics of our prob- 
lem is thereby unchanged by adding to 7i a penalty term 



T-i' = -a 



J2- J(J+1) 



with a — !■ +00 and considering 

only the 2J + 1 lowest levels. 

To study the exchange splitting between any pair of 
states, we need the following matrix element of the 
imaginary-time propagator between two coherent states: 
A = (JMn|e-^^ \JM'n') (we set h = 1). The coher- 
ent states are defined as usual by rotating a state \JM) 
from the z axis to n = (sin 6* cos sin sin 93, cos 6') , i.e., 



IJMn) 



\JM). Clearly, A is unaf- 



(2J-2j)!(./+Af)! 

{2jy.ij+M-2jy. 
equality is an exact result 



fected by the penalty term Ti.' and we can simply omit it 
for our problem. 

One can show that {JMn\ [\jjn) ^ \ J — j, M — j, n)] = 

(1 - ^-^"^JT^) ' (where the first 
and the second one an 
approximation valid for j(J — A/) ^ 2J), so that 
\JMvl) Ri Ijjn) (g) \J — j,M — j,n) to relative order 
The proof of above equality, to be detailed else- 



2J 



where, uses the fact that states of spin J can be expressed 
as completely symmetrized (over all possible permu- 
tations) tensorial products of 2J spin-1/2 states, and 
exploits the group structure of permutations. By a simi- 
lar argument (together with the fact Ti., depending only 
on the total spin J, commutes with the permutation op- 
erator), one can also prove the following exact result: A^ 

^^mJ^^ m-^\^{J-J.M+jM]e-nJM'n'). 
Combining those results, we obtain: A oc 



those 

n|(g){^JMn c^^^ IjjnO^i JMn'^ , with 

j = {M' — M)/2 (without restriction, we assume 
M' > M), J = J-j, and M = {M + M')/2. Now, 
for large values of J and small values of the ap- 
plied field along the hard axis, n and n' remain very 
close to z and — z, respectively, so that we can write 
|jj,-n) « e-|j,-j) and « e'"'|j,-j). This 

finally gives A cx J 'Du{T)e~^^^'^'^^\ where the path 

integral is for coherent states JMuj with u(0) = n 

and u(T) = n', and where the action is given, as 
usual, by 5[u(r)] = 5wz[u(t)] + Sh[u{t)]. The first 
term is the Wess-Zumino (or Berry phase) action, 
'5wz[u(t)] = iM J {1 — cosOu) difu, responsible for the 
quantum interferences |5| ; the second term is the dynam- 
ical action, 5/f[u(T)] = dr_E(u(r)), where the energy 



is £'(u) = |^(j, -j\ (g) {^JMu J n -j) (g) JMu 
In short, we have mapped our original problem onto 



JMn) and 



that of the tunnelling between the states 



JMxi'i of a fictitious spin J with biaxial anisotropy, 
which can be treated by the instanton method as in 
Ref. 0. Noting that this fictitious spin is subject to 
the effective field = — 2jK, along the easy axis, 
we immediately generate the complete set of diabolical 
points (|lal b'l. with = K; this agrees well with the 
exact result = V K'^ — D"^ for D <^ K, which is 
actually the case for Feg. By mapping the original 
tunnelling problem onto that of a fictitious spin in an 
effective Hz field, we obtain a simple interpretation of all 
the diabolical points in terms of destructive interferences 
due to the Berry phase for the effective spin J. The 
striking parity alternation discovered by Wernsdorfer 
and Sessoh Q (red vs. blue points in Fig. ^), is thus 
simply interpreted as due to J being alternately integer 
and half-integer. 

Note, that in principle, our approach is supposed to be 
valid only in the limit of large J and for small and 
Hz] it is thus a surprise to see that it essentially yields 
exact results, even for small J and/or for lar ge H x and 
Hz- This puzzle has already been noticed [loLlllj and is 
not fully understood. 
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